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ABSTRACT 

The mere location of a Beat Cepheid model in a Period Ratio vs. Period diagram (Petersen diagram) 
puts very tight constraints on its metallicity Z. The Beat Cepheid Peterson diagrams are revisited with 
linear nonadiabatic turbulent convective models, and their accuracy as a probe for stellar metallicity 
is evaluated. They are shown to be largely independent of the helium content Y, and they are also 
only weakly dependent on the mass-luminosity relation that is used in their construction. However, 
they are found to show sensitivity to the relative abundances of the elements that are lumped into 
the metallicity parameter Z. Rotation is estimated to have but a small effect on the 'pulsation 
metallicities'. A composite Petersen diagram is presented that allows one to read off upper and lower 
limits on the metallicity Z from the measured period Pq and period ratio Pig. 

Subject headings: (stars: variables:) Cepheids, stars: oscillations (including pulsations), stars: rota- 
tion, galaxies: abundances 



1. INTRODUCTION 

The idea of using the periods extracted from the 
lightcurves of Beat Cepheids (BCs) to constrain the as- 
trophysical parameters of these Cepheids is certainly not 
new. It goes back to the introduction of Period Ra- 
tio vs. Period diagrams, or Petersen (1973) diagrams 
(PDs). The usefulness of PDs comes from the fact 
that the period ratios Pio = Pi/Po turn out to be 
strongly dependent on the metallicity Z. This was sub- 
sequently exploited to show that the newly introduced 
OPAL opacities (Iglesias & Rogers 1996) largely removed 
a longstanding mass discrepancy problem {e.g., Moska- 
lik, Buchler & Marom (1992)). With the microlensing 
projects a large number of extra-galactic BCs were dis- 
covered, and the consequent utility of PDs, both ob- 
servationally, e.g., (Andrievsky et al. 1994; Beaulieu, 
Buchler & KoUath 2001), and theoretically, e.g., (Mor- 
gan & Welch 1997; Christensen-Dalsgaard & Petersen 
1995; Kovacs 2000) was quickly realized. 

The modeling of galactic evolution relies heavily on 
observational galactic metallicity distributions. The lat- 
ter have traditionally been estimated by spectroscopic 
means of a number of different objects: through the ob- 
servations of HII regions {e.g., Vi'lchez et al. (1988), Ur- 
baneja et al. (2005)), of B type supergiant stars (Mon- 
teverde et al. 1997) of Planetary Nebulae (Magrini et al. 
2004), of Wolf Rayet stars (Abbott et al. 2004), of red gi- 
ant branch photometry (Tiede et al. 2004) and of classi- 
cal Cepheids (Andrievsky et al. 1994). Techniques based 
on multiband photometry of regular Cepheid lightcurves 
(Caputo et al. 2001) and the use of the shapes of Bump 
Cepheid lightcurves (Keller & Wood 2006) have also been 
proposed. All these methods have varying degrees of 
accuracy and, in the case of M33, for example, a con- 
troversy has arisen between the 'traditional' large value 
of the [O/H] gradient ~ —0.11 dex/kpc (Vilchez et al. 
1988; Garnett et al. 1997) and a revised shallower one 
of -0.012 dex/kpc of Crockett et al. (2006). 

It is therefore extremely useful to have an additional, 
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independent and accurate method for obtaining a metal- 
licity. In this spirit Beaulieu, Buchler, Marquette, Hart- 
man & Schwarzenberg (2006) used 5 newly discovered 
BCs in M33 as metallicity tracers. From the 'pulsational' 
metallicities that they calculated and from the galacto- 
centric distances of these Cepheids they could extract a 
metallicity gradient of -0.16 dex/kpc for M33 which al- 
lowed them to weigh heavily in favor of the traditional 
value of [O/H] and against a recent downward revision. 

In this paper we explore in some detail both the power 
and the limitations of the 'pulsational' determination of 
the metallicities of BCs on the basis of linear nonadia- 
batic convective stellar models. In a subsequent paper 
we will present an extended survey of nonlinear models, 
and we show how they can be used to further narrow 
down the metallicity determination of BCs. 

2. LINEAR PULSATION MODELS FOR BEAT CEPHEIDS 

Only the envelope of the Cepheids is engaged in pul- 
sation. This allows theorists a convenient simplification: 
One needs to construct Cepheid envelopes only on which 
one imposes an inner boundary condition (v=0) at a fixed 
core radius R — R^, where in addition the luminosity is 
constant L — and is provided by the nuclear burning 
in the stellar interior, a region which is irrelevant inso- 
far as the pulsation is concerned. The envelope which 
can be taken as chemically uniform because of a prior 
evolutionary stage of extended convection, is uniquely 
specified by 3 astrophysical parameters, e.g., the stellar 
mass M , the luminosity L and the effective temperature 
Toff of the equilibrium model {i.e. the state it would be 
in if it did not pulsate), as well as the composition pa- 
rameters, i.e. X, Y and Z, that appear in the equation 
of state and in the opacities. We recall that X and Y are 
the mass fractions of hydrogen and helium, and that all 
other elements are lumped into the 'metallicity' Z. The 
three parameters, M, L, Tcs, are the most convenient 
from a computational point of view, and are the input 
to our equilibrium model builder and linear pulsation 
codes {e.g., Yecko, KoUath & Buchler (1998)). Other 
astrophysical parameters can be substituted, but at the 
expense of an iteration, e.g., the stellar radius of the 
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Fig. 1.— Loci of potential BCs as a function of Z for specified pairs. Left: (Po.Ao) = (6.188, 0.7005), (3.827, 0.7091), (1.250, 0.7250); 
right: = (6.00, 0.7150), (3.00, 0.7250), (1.250, 0.7400). Z increases by 0.001 between the successive loci. The vertical Unes represent M-L 
relations that are derived from Pietrinferni et al. (2006) for Z=0.001, 0.004, 0.008, 0.019, 0.027, from loft to right. 
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Fig. 2. — Location of the Cepheid models with simultaneously linearly unstable F and Ol modes; the slanted lines are at 
constant period, with Po ranging from 1.0 to 6.5 d in steps of 0.5 d; superposed are the Z=0.004, 0.008 and 0.019 evolutionary 
tracks of Girardi et al. (2000) for masses ranging from 3 to 6Mq. 
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Fig. 3. — Pi/Po vs. Log Pq plot. The pairs of lines delimit 
the ranges for which both F and 01 are linearly unstable. 
The metallicity increases downward in the figure from Z = 
0.001, 0.002, 0.003, 0.004, 0.008, 0.010, 0.016, 0.020, 0.026. 
For reference we have superposed the location of known Beat 
Cepheids (cf text). 

equilibrium model, R^{L, M^T^e), or one or more of the 
linear periods, Pk{L, AI,Tcs), where k = 0,1,... label 
the radial pulsation modes. Note that the three basic 
parameters are also the natural quantities L{t), M{t), 
and Teg(t) whose temporal behavior, i.e. the evolution- 
ary tracks, are given by stellar evolution calculations. 

Our stellar equilibrium models are constructed and 
their linear (nonadiabatic) periods Pk and growth rates 
rjk are computed with the Florida pulsation code. Tur- 
bulent convection is approximated with time-dependent 
mixing length (described e.g., in KoUath et al. (2002, 
1998)) e.g., with the values of the a parameters chosen 
to be {ad= 2.177, aa= 0.4, as= 0.433, a„= 0.12, at = 
0.001, ar= 0.4, ap= 0, ax^ 1.5}). OPAL (Iglesias & 
Rogers 1996) and Alexander & Ferguson (1994) opaci- 
ties are used. 

For our purposes we can use linear (nonadiabatic) pe- 
riods. Linear periods are sufficiently accurate because 
they differ from the nonlinear ones at most in the fourth 
decimal figure {e.g., Antonello & Aikawa (1998), Szabo 
& Buchler (2007)). It is in the modal selection problem 
where nonlinearity plays an important role in limiting 
the region where beat pulsations can occur (sec KoUath 
& Buchler (2001), KoUath et al. (2002)). We wiU ad- 
dress this issue in Szabo & Buchler (2007) which deals 
with nonlinear BC model pulsations. 

During the galactic evolution metal enrichment occurs 
through helium and hydrogen burning. Because the he- 
lium content Y thus correlates with Z , we have assumed 
a relation X = X{Z), obtained by fitting a parabola in 
Log Y vs. Log Z through the 3 points {Y , Z)= (0.280, 
0.020), (0.276, 0.008) and (0.270, 0.004). With this re- 
lation we reduce the 3 composition parameters {X{Z), 
Y{Z)=l-X{Z)-Z, Z) to a single one, Z. 



3. RESULTS 
3.1. Log M - Log L plane 

In their calculation of the metallicities of the BCs in the 
M33 galaxy (Beaulieu et al. 2006) used the 2 observed 
periods Po{L, M, T^g, X, Z) and Pi{L, M, T^g, X, Z), and 
the Toff as the 3 parameters that determine the stellar 
model. This then produced, for a given composition, a 
one parameter family of models, with T^g as the param- 
eter. They then further limited the models to those that 
were linearly unstable simultaneously in the fundamental 
F and first overtone Ol modes. This yielded, for each 
Z, loci of potential BCs in a Log L vs. Log M diagram 
(potential because this double instability is only a neces- 
sary condition for BCs to exist). 

In Fig. 1, to illustrate the procedure, we similarly dis- 
play the Z dependent set of loci of 6 hypothetical BCs, 
with the {Pq,Piq) values that are indicted in the figure. 
For each sequence (in Tcs, with fixed Pq and Pio) we 
have plotted the loci as a function of Z. For the top 2 
sequences on the left, Z ranges from 0.003 to 0.017, and 
for the other four it ranges from 0.001 to 0.020. 

Superposed on the loci we display the M-L relations 
derived from Pietrinferni et al. (2006), for Z= 0.001, 
0.004, 0.008, 0.019 and 0.027, from left to right. The 
acceptable BC candidate models are those for which the 
locus of a given Z intersects the M-L curve for the same 
Z. As noted in Beaulieu et al. (2006), as Z increases 
the loci move toward the northeast, at first very slowly, 
then they speed up as they cross the region of the M-L 
curves to slow down again when they arrive on the right. 
It is this rapid crossing of the AI-L region that allows 
one to narrow down the Z of the BCs. We recall that 
the loci depend only on the period Pq and period ratio 
Pio, and pulsation theory. Stellar evolution comes into 
play through the M-L relation which determines what 
are the acceptable metallicity values for the given Pq and 
PiQ. It is apparent from Fig. 1 of Beaulieu et al. (2006) 
that these values are relatively insensitive to the chosen 
M-L relation (see also §3.4 below). 

While this procedure allows us to determine the metal- 
licity of BCs, it is computer intensive as it necessitates 
the computation of the loci for a range of Z, and this for 
each BC. Furthermore for some parameter values that 
fall far away from those of normal Cepheids the models 
do not always converge. 

3.2. Log Pq - Piq plane (Petersen Diagrams) 

We find it useful to present a different, more practical 
procedure. This consists of constructing Petersen dia- 
grams (Pio , Pq] Z) in which we delimit the ranges of ac- 
ceptable BC models for a finely spaced set of Z values. 
For a given BC with its observationally derived Pq and 
Pio one can then readily determine into which of these 
ranges it falls, interpolate, and thus impose upper and 
lower limits on Z. 

Our PDs are constructed as follows. We specify the 
F period Po{L, M,Tes, Z), and the M-L relation L = 
L{M,Z). These 2 constraints allow us to compute a 
sequence of models of a specified Z with different values 
of Tcff and consequently different Pio(Af, L, Tcff, ■^). Of 
these models we then only retain those that arc linearly 
unstable both in the F and the 01 mode. Our set of 
Pq ranges from 0.75 d to 6.5 d. Our Z values range from 
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Fig. 4. — Log Po vs. Pi/Po {Po on top axis). The lines delimit the ranges for which both F and 01 are linearly unstable: 
Upper limit on the left and lower limit on the right. The metallicity increases downward in the figure from Z = 0.001 (top line) 
to 0.010 in steps of 0.001, and from 0.010 to 0.026 (bottom line) in steps of 0.002. (A color figure with AZ =0.001 spacing 
everywhere is available with the electronic version.) For reference we have superposed the location of known Beat Cepheids. 
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Fig. 5. — Period ratio vs. Log Tcb for models with simultaneously linearly unstable F and 01 modes, for 3 different compositions. 
The lines correspond to models with constant periods, running from 0.5 to 6.5 days. 
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Fig. 6. — Period ratio Pio vs. mass at constant period for 
BCs candidates with the various values of Z from the top 
down as indicated in the figure. The periods range from 1.0 to 
6.5 days, with the larger masses going with the larger periods. 

0.001 to 0.026 in steps of AZ = 0.001. 

We recall that there exists a slightly slanted, vertical 
strip in the HR diagram, called the instability strip (IS) 
in which the stars become unstable to pulsation and are 
called Cepheids. In Figure 2 we display our doubly un- 
stable models in an HR diagram for selected metallicities 
(Z=0.004, 0.008 and 0.019). The lines run at an angle 
because the sequences (in Toff) are computed at constant 
period Pq and constant Piq. The full IS is much wider in 
Toff as its extends leftward and rightward into the regimes 
where onij the 01 and F modes, respectively, are linearly 
unstable. Note also that not the whole doubly unstable 
region will give rise to BCs: first, simultaneous instability 
in the two modes is only a necessary condition for beat 
behavior, and second, the tracks that actual Cepheids 
follow may indeed not penetrate everywhere into this re- 
gion. In these calculations we have made use of a M-L 
relation derived from Girardi et al. (2000) shown in Ta- 
ble. 1 (we will return to a discussion of M-L relations in 
§3.4). 

In Figure 2 we superpose the tracks of Girardi et 
al. (2000) for Z=0.004, 0.008 and 0.019, respectively. 
Among the many tracks in the literature we have cho- 
sen these because they cover the whole mass range and 
Z range of interest to us. Furthermore we have found 
that the pulsational properties of the Cepheid models 
along these evolutionary tracks, as calculated with our 
pulsational code, are in good agreement with the reso- 
nance constraints imposed by the OGLE-2 data of the 
Small and Large Magellanic Clouds (Buchlcr, KoUath & 
Beaulieu 2004). 

The key results of our survey are displayed as a PD in 
Fig. 3. The pairs of lines delineate for each Z the range 
of Pio values for which the F and 01 modes are simulta- 
neously linearly unstable. The mctallicity is increasing 
downward for selected values of Z = 0.001, 0.002, 0.003, 
0.004, 0.008, 0.010, 0.016, 0.020, 0.026. The lines be- 
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Fig. 7. — Pi/Po vs. Log Po plot. Effect of changing the 
helium content Y {=1-X-Z) at fixed Z. The pairs of lines 
delimit the ranges for which both F and 01 are linearly un- 
stable. 

come crowded for the higher Z indicative a spacing that 
is closer to Log Z than to Z. 

Because of overcrowdedness we could only display very 
few values of Z in Fig. 3 which we have therefore split 
into 2 subfigures Fig. 4, showing respectively the lower 
and upper limits of Pio for a finer set of Z values. The Z 
are increasing downward from Z=0.001 to 0.010 insteps 
of 0.001, and then from Z = 0.010 to 0.026 in steps of 
0.002. 

For reference we have superposed the location of 
known BCs in the Galaxy (x's), in M33 (Beaulieu et al. 
2006) (hexagons), in the LMC (triangles) (Alcock et al. 
1995; Soszyhski et al. 2000) and in the SMC (crosses) 
(Beauheu et al. 1997; Udalski et al. 1999). 

Fig. 3 and, more conveniently. Fig. 4 can now be used 
to infer upper and lower bounds on the metallicity of a 
given BC star from its coordinates, (LogPg, Pio)- We 
identify the range(s) into which a star falls which gives a 
Zrnin and Z„iax- If it falls between the calculated ranges, 
as happens for short periods, we get a Zmin and Z^ax 
from the adjacent ranges. One can refine Zmin and Zmax 
by interpolation in the table from which the figures are 
constructed. 

Consider, for example, the M33 BC (hexagon) with the 
lowest PiQ. The left panel of Fig. 4 shows that the star 
falls between the Z=0.010 and 0.011 lines, and from the 
right we find 0.012< Z <0.013. An interpolation gives 
0.0106 < Z < 0.0124 with an average of Z=0.0115. 

The 'Beat Cepheid metallicities' that we infer for the 
Galactic BCs, Z=0.0074 ~ 0.0182, with an average of 
Z=0.0118 are on the low side, but for the MC and the 
M33 BCs they fall in the generally accepted ballpark, 
0.001 - 0.007 with average Z = 0.0037 for the SMC, 
0.0035 - 0.012, with average = 0.0062 for LMC, and 
0.0075 - 0.0124, with average = 0.0075 for M33. We 
will return to a discussion of a possible remedy for this 
discrepancy. 
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Fig. 8.— M-L relations: Z =0.001, 0.004, 0.008, 0.019, from 
the top down. 

In the upper left corner of Fig. 3 , the ranges for the 
plotted Z values do not overlap, meaning that the posi- 
tion in the PD very narrowly determines the allowed Z 
for a given Pg ^-nd Piq. In the upper right, in contrast, 
there is more overlap and the allowed range of Z is a 
little broader, albeit still in the ±0.001 range for Z. In 
general toward the bottom there is less overlap as well. 

Fig. 4 provides us with no information about how the 
temperature or mass vary between the upper and lower 
loci at fixed Pq. In Fig. 5 we display the behavior of the 
period ratio Pio as a function of T^s, each sequence at 
fixed period Pq. We see that for high Z (on the right) 
the upper boundary, i.e. the maximum Pio, of the doubly 
unstable models is at a higher Toff for the short periods, 
but the reverse is true for the longer periods. For the low 
Z (left) the behavior of Pio with Toff is not even mono- 
tone. Referring back to Fig. 3 it is therefore apparent 
that the behavior of Toff between the boundaries of the 
doubly unstable models (vertically, at fixed Pq) can be 
increasing, decreasing or even be nonmonotone. Fig. 6 
indicates that similar nonmonotone behavior occurs for 
Pio as a function of mass. 

In the following we examine what physical factors can 
influence the PD, and we will examine in turn the effects 
of the helium mass fraction Y of the M-L relation, and 
of stellar rotation. 

3.3. Effect of the Helium Content Y 

Because we have chosen the helium content relation 
Y{Z) in a somewhat cavalier way, we want to explore 
the sensitivity of the ranges of BC models in the PD to 
Y . In Fig. 7 we show the ranges for 3 selected values 
of Z = 0.004, 0.008 and 0.20. The figure also displays 
the associated X and Y values. It is very clear from this 
figure that the sensitivity to even large changes of Y is 
minimal. The period ratios are essentially only sensitive 
to Z (because opacity mostly is mostly sensitive to Z). 

3.4. Effects of the M-L relation 
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Fig. 9.— P\/Po vs. Log Pq plot. Effect of the M-L relation: 
Girardi et al. (2000) (black), Pietrinferni et al. (2006) (red) 
and Bono et al. (blue). The pairs of lines delimit the ranges 
for which F and Ol are simultaneously linearly unstable. 

It is well known that a Cepheid star of given mass may 
cross the IS three times, as seen e.g., in Fig. 2 for the Gi- 
rardi et al. (2000) tracks. The first crossing occurs very 
fast, and, statistically, very few stars would be caught in 
that state of evolution. Most Cepheids are expected to 
be found on the slow leftward (blueward) loop, and some 
on the faster rightward return. 

Clearly the tracks of Fig. 2 show some serious and 
well known deficiencies as far as observation is concerned 
(which are shared by all other published evolutionary 
tracks): the blueward 'Cepheid' loops peter out as Z is 
increased, and by the time one reaches Galactic values of 
Z, there are no more loops at all that penetrate the BC 
Cepheid region, in strong discordance with observation. 
While some decrease of the lengths of the horizontal loops 
with M is acceptable because there is a lower observed 
cutoff in the Cepheid periods, it is patently excessive. 
A good discussion of the problems associated with the 
horizontal Cepheid loops can be found in Cordier (2000) 
and Cordier et al. (2002). 

Aside from these problems, the very use of an AL-L 
relation is an approximation, first, because we limit our- 
selves to the second crossing, and second, because even 
then the stellar tracks are not horizontal in the HR dia- 
gram (L-Toff plane). In addition, referring back to Fig. 2 
one realizes that extracting an M-L relation from the 
evolutionary tracks requires some guesswork as to where 
some of the aborted horizontal loops really should be. 
Another problem is that in essentially all the evolution- 
ary calculations the horizontal loops exhibit a nonmono- 
tone behavior with Z, some of which may be of numerical 
origin. There is also broad evidence from the LMC and 
SMC Cepheids that the computed evolutionary tracks 
may be undcrluminous for a given mass to account for 
their pulsational properties {e.g., Beaulieu, Buchlcr & 
Kollath (2001)), and a smilar conclusion was reached for 
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TABLE 1 
M-L RELATIONS 
LogL 

■ + b2LogM + b^LogZ + biLog'^M + b^Log'^Z + bgLogZLogM 





bi 


62 63 


b4 


b5 


fee 


Gir 


-3.2822 


8.0864 -1.7988 


-1.2764 


-0.1771 


0.9588 


Piet 


-2.0937 


5.5361 -1.5061 


0.0306 


-0.1421 


0.7112 



the bump Cepheids by (Keller & Wood 2002; Bono et 
al. 2002). However, it is not our purpose here to present 
a critical review of the various evolutionary calculations 
and the M-L relations one can derive from them, espe- 
cially, as it will turn out that our application to BCs 
fortunately has only a moderate sensitivity to the M- 
L relation. This was already noted in Beaulieu et al. 
(2006) and in Fig. 1 above. 

With these caveats, using the blue tips of the Gi- 
rardi et al. (2000) tracks we have made fits to Log L 
that are bi-quadratic in Log M and Log Z in the range 
2.5 < M < 6.5 and 0.004 < Z < 0.02. They have the 
largest uncertainty for large Z in general, and for small 
M (because of the lack of horizontal loops in the evolu- 
tionary tracks) . Table 1 presents the M-L relations that 
we have extracted from Girardi et al. (2000) and from 
Pietrinferni et al. (2006). Both these evolutionary cal- 
culations give a full grid of M values and Z values in the 
range that we are interested in, and they use a nonzero 
overshooting parameter. Alibert et al. (1999) give an 
M-L only for some values of Z and we have therefore not 
used them here. Cordier et al. (2002) find that the evo- 
lutionary calculations without overshooting does not give 
good agreement with observational constraints. In addi- 
tion, we find that for Cepheids with periods as low as 15 d 
with Bono's M-L relation one would need to go to very 
high masses of 8.5 - 9 Mq as opposed to 7.0 - 7.5 Mq for 
the other M-L relations. However, just for comparison 
we also consider the M-L of Bono et al. (2000) derived 
from tracks without any overshooting. M-L fitting for- 
mulae are presented in their paper. 

Fig. 8 displays the M-L relations of Girardi et al. 
(2000); Pietrinferni et al. (2006); Bono et al. (2000) 
for Z = 0.001, 0.004, 0.008 and 0.019. One sees that 
Bono is much flatter than the other two. 

The period ratios Pio obtained with the three M-L re- 
lations are compared Fig. 9. The metallicities range from 
Z=0.001 (top), 0.004, 0.008 to 0.020 (bottom). Despite 
the large difference in the 3 M-L relations, one finds, 
perhaps unexpectedly, that the PDs for the 3 relations 
agree quite well for periods down to Pq = 1.5 d for all 
Z . The discrepancies between the M-L relations show 
up for short periods only, and especially for Bono. 

3.5. Effects of Rotation 

The effect of rotation on PDs has been considered in 
the past, generally with the help of the crude, but rea- 
sonable simplifying assumption of a 'spherical centrifugal 
force', F = uf^r (Stothers, R. (1974); Kovacs & Buchler 
(1990)) that is added in the model builder and in the lin- 
ear stability analysis. Note that this subterfuge is likely 
to exaggerate the effect of rotation. 

The results of our calculations with pseudo-rotation 
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Fig. 10.— Pi/Po vs. LogPo plot. Effect of the rotation. The 
pairs of lines delimit the ranges for which both F and Ol are 
linearly unstable. The curves are for Vrot = 0, 10 and 20 km/s. 
Note that the effect of rotation is most likely overestimated 
by the treatment as a pseudo-rotation. 
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Fig. 11. — Pi/Po vs. Log Po plot. Effect of changing the 
chemical makeup of Z: heavier elements reduced by a factor 
of 0.75 (blue), which shifts the curves upward to higher Pio. 

are displayed in Fig. 10 for surface rotation velocities of 
Vrot — 0, lOkm/s and 20 km/s. The PD indicates only 
small shifts even with rotations as high as 20 km/s. A 
comparison with Fig. 3 shows that for the low Z mod- 
els, with their high Pio values, the neglect of rotation 
amounts up to 20 km/s generally leads to an error in Z 
of less than 0.001, but up to 0.002 for the longer periods 
and smaller Pio (larger Z). 
Recently a flag has been raised by Suarez, Garido & 
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Goupil (2006) who found that the PDs of delta Scuti 
models exhibit a large sensitivity to the assumed rota- 
tion rate, and that even a modest rotation rate can cause 
an increase in the period ratio which would then be mis- 
interpreted as a lower metallicity than the star has. In 
contrast, we have seen that for Cepheid models rotation 
has only a small effect on PDs, all the more so that the 
observed Cepheid rotation rates (u, ot sin i) are typically 
less than 20km/s (Nardetto et al. 2006) and thus much 
smaller than those of delta Scuti stars. 

We conclude that neither the helium mass fraction Y , 
nor the stellar rotation, nor the M-L relation can appre- 
ciably affect the PD. Remaining uncertainties are asso- 
ciated with the relative elemental concentrations of the 
'metals' in Z . 

3.6. Effect of Elemental Mixture that is Lumped into Z 

So far we have used the OPAL (Iglesias & Rogers 1996) 
opacities with the standard solar elemental mixture of 
Grevesse & Noels (1993). One needs to ask the question 
of whether the PDs are sensitive to Z as a whole, or 
to the relative fraction of the heavier metals, i.e. the Fe 
group elements, which provide the bulk of the opacity. 
We therefore have made several tests. 



In Figs. 11 and 12 we display PDs obtained with OPAL 
opacities in which the number densities of the dominant 
heavy elements (Fe, Ca, Ar, S, Si and Mg) are arbitrar- 
ily lowered by factors of 0.75 (before renormalization to 
unity) in a 'Type 1 OPAL Table' with fixed X, y and Z. 
This can be compared to the PD of Fig. 3 which was for 
a solar mix (GN93 in OPAL). One notes that, in first ap- 
proximation, the lines are just shifted, but there is also a 
deformation. For example, the old Z=0.022 line in Fig. 3 
is now the Z=0.026 hne in Fig. 12, and the Z=0.005 is 
close to the 0.006 line. 

As far as the Galactic BCs are concerned we find a 
range of Z ^ 0.0089 - 0.0221, with an average of Z= 
0.0141, better than the average of 0.0118 that we find 
with a solar mix, but still a little low compared to the 
generally quoted values. For the SMC we have Z = 
0.0015 - 0.0100 with average = 0.0047, for the LMC 
0.0047 - 0.0148, average = 0.0076, and for M33 0.0013 - 
0.0146, average = 0.009. 

We are led to conclude that there is a discrepancy be- 
tween the stellar models that use opacities with the GN95 
solar mix and the metallicities that are quoted in the 
literature. Better agreement can be achieved with a re- 
duction in the relative number densities of the Fe group 
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elements. 

We have also examined other changes. The Grevcsse 
G91 mix (in OPAL) has a reduction of about 10% in the 
heavier elements. Not astonishingly, the results fall in be- 
tween those of the 25% reduction and GN95. Similarly, 
increasing the heavy elements (Fe, Ca, Ar, S, Si and Mg) 
by a factor of 1.1 produces a comparable change, but in 
the opposite direction, as expected. The alpha-enhanced 
opacities, labeled Charboyer in the OPAL opacity li- 
brary, which have a chemical makeup that is quite dif- 
ferent from solar produce large changes in the PD. Space 
docs not permit us to look with more detail into the ef- 
fects of the chemical makeup in order to see if we can 
isolate what effects the various subgroups of elements in 
the mix have on the PDs. This will be done in a subse- 
quent paper. 

We conchide that the chemical makeup does play a 
non-negligible role, mostly for the short periods. De- 
creasing the relative Fe content at fixed Z is seen to be 
equivalent to decreasing Z at fixed elemental mixture in 
first approximation. This confirms what one may have 
expected, namely that the PD are predominantly sensi- 
tive to the strength of the opacity which is provided by 
the Fe group elements. 

4. CONCLUSIONS 

With the help of our convectivc pulsation code we have 
constructed Petersen diagrams (Pio vs. Log Po) for Beat 
Cepheids. In these PDs we have delineated the region 
of simultaneous linear instability in the F and the 01 
modes for a range of relevant values of metallicity Z. 
Our results are presented in Fig. 4 such a way that from 
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the measured period Pq and period ratio Pio one can 
determine upper and lower limits on the metallicity Z of 
the corresponding Beat Cepheid. 

The results are found to be rather insensitive (1) to 
the assumed helium content Y(Z), (2) to the exact form 
of the mass-luminosity relation, and (3) to the presence 
of stellar rotation. By comparing PDs with solar abun- 
dances to those in which the heavier elements (in partic- 
ular Fe) are reduced or increased relative to the lighter 
elements we conclude that the PDs are sensitive to the 
chemical makeup of Z. 

Interestingly, we find that with a reduction in the num- 
ber densities of heavier elements at fixed Z the 'pulsation 
metallicities' for the Galactic Beat Cepheids are in better 
agreement with the values that are in the literature than 
with the solar elemental mix of Grevesse-Noels (1993). 
We plan to examine the effect of the relative abundances 
in greater detail in a companion paper, especially in light 
of recent suggested revisions of the abundances (Asplund, 
Grevesse & Sauval 2005). 
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